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ABSTRACT 

We study extremal black hole solutions of D = 5 Gauss-Bonnet gravity cou- 
pled to a system of gauge and scalar fields. As in Einstein gravity, we find that 
the values of the scalar fields on the horizon must extremize a certain effective 
potential that depends on the black hole charges. If the matrix of second deriva- 
tives of the effective potential at this extremum has positive eigenvalues, we give 
evidence, based on a near horizon perturbative expansion, that the attractor 
mechanism continues to hold in this general class of theories. We numerically 
construct solutions to a particular simple single scalar field model that display 
the attractor mechanism over a wide range of asymptotic values for the scalar 
field. We also numerically construct non-extremal solutions and show that the 
attractor mechanism fails to hold away from extremality. 



manber@physics . umass . edu 
2 kastor@physics . umass . edu 



1 Introduction 



The attractor mechanism was originally discovered for BPS extremal black holes in M = 2 
supergravity theories [1, 2, 3]. Under the attractor mechanism, the values of moduli scalars at 
a BPS extremal black hole horizon are fixed, independently of their values at infinity, in terms 
of the electric and magnetic charges carried by the black hole. More recently, begining with 
the work of references [4] and [5], the attractor mechanism has been investigated for extremal 
black holes in non-supersymmetric theories, as well as for non-BPS extremal solutions in 
J\f > 1 supersymmetric theories. 

The methods employed in references [4] and [5] are quite different and offer complemen- 
tary insights into the physics of the attractor mechanism. Reference [4] focuses on the near 
horizon limit, which is assumed to have geometry AdS2 x S D ~ 2 . An "entropy function" is 
defined by taking the Legendre transform, with respect to the electric charges, of the integral 
of the Lagrangian density over the S D ~ 2 . The constant values of the moduli fields, which 
solve the equations of motion in the near horizon region, can be shown to be those which 
extremize this entropy function. These are the attractor values of the moduli fields. Further, 
the black hole entropy is given by the value of the entropy function at its extremum. The 
approach of reference [4] is quite general and by design includes the possibility of higher 
derivative gravitational interactions. For instance, equality of the extremum value of the 
entropy function with the black hole entropy is demonstrated by showing that it reproduces 
Wald's formula [6, 7, 8, 9] which holds in higher derivative gravity theories. 

In contrast, reference [5] focuses only on Einstein gravity, but obtains very explicit results 
that illustrate the attractor mechanism in action. The authors use a combination of analytic 
and numerical techniques to follow the radial evolution of the moduli scalars from their 
attractor values, at an extremal black hole horizon, out to independent values at infinity. 
Although this approach lacks the generality of reference [4], one sees the operation of the 
attractor mechanism in a vivid way. 

In this paper, we will follow the approach of reference [5] and study the attractor mech- 
anism with a simple higher derivative gravitational interaction, the Gauss-Bonnet term, 
added to the action. We focus on D = 5, the smallest dimension in which the Gauss-Bonnet 
interaction is non-trivial. This theory serves as the first example, beyond Einstein gravity, 
of a Lovelock gravity theory [10]. Lovelock gravity theories share a number of important 
properties with Einstein gravity [10] and have been studied in many contexts over the years. 
In particular, vacuum and electrovac black hole solutions have been well studied, begining 
with the work of [11, 12, 13]. 

One knows from the general results of [4] that the moduli scalars in Gauss-Bonnet gravity 
must take values at an extremal black hole horizon that extremize the entropy function for 
this theory 3 . Our results establish, at least within the particular gauge and scalar field system 
we study, that these near horizon attractor values are actually obtained at the horizon in 
asymptotically flat, extremal black hole solutions with a range of different values for the 
scalars at infinity 4 . 

3 Thc entropy function for Lovelock gravity theories, which includes Gauss-Bonnet gravity, is calculated 
explicitly in reference [14] (see also [15]). 

4 Non-supersymmetric, asymptotically flat extremal black hole solutions in a higher derivative gravity the- 
ory have also been studied in reference [16]. The theory considered in this case is a general R 2 gravity theory 
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In section (2) we describe the D = 5 Gauss-Bonnet gravity theory coupled to a system 
of gauge and scalar fields that we will be studying. In section (3) we present a set of simple 
analytic black hole solutions in which the scalar fields take constant values throughout the 
spacetime. The possible constant values are, as in reference [5], the extremal points of a 
certain effective potential function that depends on the charges carried by the black hole 5 . In 
section (4) we numerically construct extremal solutions in a certain single scalar field model, 
in which the scalar fields vary between fixed attractor values as the degenerate horizon and 
independent values at infinity. These attractor values are again the extremal points of the 
effective potential, with the provision that the eigenvalues for small fluctuations about the 
extremal point must all be positive for the attractor mechanism to hold. We evaluate the 
ADM mass of these spacetimes and show that it is minimized by the extremal solutions of 
section (3) in which the scalars are constant at their attractor values. Finally, in section 
(5) we numerically construct non-extremal black hole solutions with non-constant values for 
the scalar field. We show that, as in the case of Einstein gravity [5, 18], non-extremal black 
holes do not exhibit the attractor mechanism, i.e. the value of the scalar field at the horizon 
depends on its value at infinity. 
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2 Gauss-Bonnet gravity 

We consider Gauss-Bonnet gravity in five dimensions coupled to N Abelian gauge fields A ^ 
with a = 1, . . ., JV and n massless scalar moduli fields <pi, with 1 = 1, ...,n. The moduli 
scalars have vanishing potential, but couple to the gauge field kinetic terms through the 
matrix function f a b(<f>)- The action is then given by 

S = ± [ dx 5 V=9 \R + aC^ -2d,<j )I dy I - f ab {ct ) j)F; v F b ^ 

K J 

with the Gauss-Bonnet term in the Lagrangian given by 

C < GB ) = R 2 - AR 1& R? S + R l5 ^R? 5Xc , (2) 

The coupling constant a has dimensions (length) 2 . The equations of motion and further 
details of this theory are given in the Appendix. 

We are interested in extremal black hole solutions of this theory. We will therefore assume 
a static, spherically symmetric form for the metric 

dr^ 

a z [r) 

where dVt 2 = r )^dx l dx^ with i,j = 1,2,3 represents the round metric on the unit 3-sphere 
in a general set of coordinates. We will restrict our attention to black holes that carry only 
electric charges for the gauge fields A®. The equations of motion for the gauge fields may 
then be solved by taking the field strengths to be of the form 

F a = ^^dt A dr. (4) 



in D = 4 with a moduli dependent coupling. Black hole solutions manifesting the attractor phenomenon are 
implicitly constructed via a series expansion method. 

5 See also reference [17] for similar results in D = 5 Maxwell-Einstein supergravity. 
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Here, the constants Q a are the electric charges and the field dependent tensor f ab (4>) is the 
inverse of the tensor coupling f a b{4>) that appears in the Lagrangian. With this form for the 
field strengths, the stress-energy tensor for the gauge fields can be written in terms of the 
effective potential 

V eS (4>) = f cd (<j>)Q c Q d . (5) 

Furthermore, the effective potential acts as a potential in the equations of motion for the 
moduli scalars, which are given by 

dr^a^dAr)-^^ (6) 

It will be important to note that these equations of motion may be solved by constant values 
4>j of the scalar fields, if these values represent a critical point of the effective potential, i.e. 
they satisfy 

(9/Veff)U = 0. (7) 

As in reference [5], we will see below that the critical points 0/ of the effective potential 
represent possible attractor values for the moduli scalars, provided that the matrix of second 
derivatives of the effective potential 

M Z j = {djdjVrf)]} (8) 

has positive eigenvalues. 

Given the form of the ansatz for the metric (3) and the gauge fields (4), the gravitational 
equations of motion are given by 

Srr = (drh) (drf) - , Q u = « 4 (fi^) (d^) + , 

<3ij = (-6 V(fl^)(fl^') + ^j^j 7 « , (9) 

where the tensor = G^ u + a Gj£ B ^ combines the Einstein tensor G^ u and its Gauss- 
Bonnet counterpart Gj^ B ^ which is given by equation (43) in the Appendix, where we also 
give expressions for the nonzero components of G^ and Gj£j B ^ in the spherically symmetric 
ansatz. 

One particular combination of the gravitational field equations will be especially impor- 
tant in the analysis below. Given the gravitational field equations in (9) and the explicit 
expressions for G^ u and Gj£j B ^ in the Appendix, one can compute the quantity g'" r Q rr — g lt Qtt 
in two ways. Setting these equal then gives the equation 



-3 a - 



b 2 + 4a(l-aV 2 )] b" = 2 a 2 (d r J )0 r 7 , (10) 



which implies that b" = if the scalar fields are constant. Note that in the absence of 
the Gauss-Bonnet term, i.e. with a = 0, equation (10) implies that b" < 0, which is an 
important ingredient in the 'c-theorem' of reference [19]. 
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3 Double-extreme solutions 



Extremal black hole solutions in which the scalar fields take constant values are sometimes 
called double-extreme solutions. In this section, we investigate such double-extreme solutions 
in the Gauss-Bonnet gravity theory of section (2). As noted above, the constant values 0/ for 
the scalar fields must be critical points of the effective potential, in order that equation (6) 
be satisfied. A second observation is that, with constant values of the scalar fields, equation 
(10) implies that b" = 0. 



Robinson-Bertotti solutions 



We first consider solutions of the Robinson-Bertotti form. The metric is taken to be AdS 2 x 
S* 3 , which we write as 

ds 2 = -^-dt 2 + ^dx 2 + b 2 H dnj. (ii) 

R l x z 

where R and b H are constants. This will be the near horizon form of the metric, both for the 
double-extreme black hole solutions which we study in this section, and the solutions with 
non-constant scalars which we subsequently study via numerical methods. 

First note that since the metric function b(r) in the Robinson-Bertotti metric is constant, 
we have b" = 0. Equation (10) then implies that the scalar fields 0/ must be constant. The 
constant values 0/ of the scalar fields must in turn satisfy equation (7). Let V be the constant 
value of the effective potential throughout the spacetime, V = V e g((f)). One can then show 
that the remaining field equations imply that the S* 3 and AdS 2 radii are given according to 

b% = ±V, R 2 = \b 2 H + a (12) 

This differs from the Einstein case only in the contribution of the Gauss-Bonnet coupling 
constant a to the AdS 2 radius. 

Extremal black hole solutions 

We now consider double extreme black hole solutions 6 . In order to satisfy b" = with 
asymptotically flat boundary conditions, we may set b(r) = r without any loss of generality. 
Plugging this into the second equation in (9) yields an equation for a(r) 



a' 



3ar 2 + 12 aa(l- a 2 )] + 3r a 2 = 3r - ^ . (13) 



We can integrate this equation to obtain the general solution 



y 2 

a 2 = ! + -- + 



r 2 V 1 / . „ , „ V 



4a \\ Aa) 2ar 2 \ 3 



r 4 - 2Mr 2 + -), (14) 



6 Because the scalar fields are constant, the black hole solutions we find here are equivalent to those of 
Gauss-Bonnet gravity coupled to a single U(l) gauge field given in references [20, 21], with the effective 
electric charge Q 2 oc V. 
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where M is a constant of integration. One can check that the remaining gravitational 
equations of motion are also with a(r) having this form. 

This general solution describes both non-extremal and extremal black holes, as well as 
naked singularities. Moreover, it is well known that Gauss-Bonnet gravity generally has two 
distinct constant curvature solutions. If the cosmological constant vanishes, as it does in the 
Lagrangian (1), then flat space is always one of these solutions. For Gauss-Bonnet coupling 
a > (a < 0) the second constant curvature solution has negative (positive) curvature. We 
need to sort through these various possibilities to identify the asymptotically flat, extremal 
black hole solutions in (14). 

Taking the large r limit, we see that for a > the + branch of (14) is asymptotically 
anti-deSitter, while with a < the — branch is asymptotically deSitter, in accordance with 
the remarks above. For asymptotically flat solutions, we must then take the — (+) branch 
of (14) for a > (a < 0). Taking the large r limit in these cases leads to 



a 2 ^l-2'—^ + — + .... (15) 



M + a V 



r 



3r 4 



Given the normalization of the Einstein term in the Lagrangian (1), the ADM mass can be 
read off from the expansion g tt ~ — 1 + Madm/Q^ 2 ^ + • • •• Comparing with equation (15) 
we see that the ADM mass of our solutions is given by Madm = 12n 2 (M + a). 

We must now analyze the horizon structure of the asymptotically flat solutions. For 
Gauss-Bonnet coupling a > 0, taking the + branch in (14), we see that a 2 = when 
the second term under the square root vanishes. Horizons then occur at the roots of the 
polynomial r 4 — 2Mr 2 + V/3. These roots are given by 

r± = M ± \J M 2 - V/3 (16) 



and we see in turn that the solutions represent black holes for M > J V/3, with the ex- 



tremality condition being M = yV/3. Keeping track of the signs carefully, one finds that 
the horizon radii are the same for a < with one proviso. In order for the metric function 
a 2 to vanish for a negative, the quantity 1 + r 2 /4a must be negative at the horizon. For the 
extremal case of interest to us, this imposes a lower bound on the Gauss-Bonnet coupling 
constant 



1 

a 



> -tV^A ( 17 ) 



4 

Only for a satisfying this bound do we find double extremal black hole solutions. 7 



Near horizon limit 

Let us now take the near horizon limit of the double extreme solutions and check that it 
coincides with the Robinson-Bertotti solutions found at the beginning of the section. In 
the extremal limit we find that the outer horizon radius r + in equation (16) is given by 
r 4 = M 2 — V/3. In terms of the parameter bn of the Robinson-Bertotti solutions, we then 

7 Extremal black hole solutions in Gauss-Bonnet gravity are also discussed in references [22, 23]. 
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have r + = bn- The near horizon limit of the metric function a 2 in (14) in this case is found 
to be 

Setting x — r — bn and i? 2 = a + b 2 H /4:, we see that the near horizon limit of the double 
extreme black hole solutions indeed coincides with the Robinson-Bertotti metric (11). 



4 Attractor mechanism 

In the last section, we found extremal black hole solutions with constant scalars in our 
theory. As in Einstein gravity [5], the constant values taken by the scalars <pi must represent 
a critical point of the effective potential V e s. In this section we will construct extremal black 
hole solutions with non-constant scalars. These will have the same near horizon limit as the 
solutions of section (3). However, we will see the asymptotic values of the scalar fields may 
be varied freely. The existence of these solutions in Gauss-Bonnet gravity establishes the 
operation of the attractor mechanism in this theory. 



A specific model 

In order to find extremal black hole solutions with non-constant scalars, we must further 
specify the theory by choosing definite numbers of gauge and scalar fields and the form of the 
couplings between them. Following reference [5], we consider a simple example that consists 

of one scalar field coupled to two U(l) gauge fields A® with a = 1,2. The couplings of the 
scalar field to the gauge fields is taken to be 

f ab (<P) = e- a ^5 ab . (19) 

It is then straightforward to compute the effective potential for this model, which is given 
by 

VM) = e a ^Ql + e a ^Q 2 2 . (20) 

In order that a critical point of V c s should exist, the constants a± and «2 must have opposite 
signs. The critical value of the scalar field is then given specified by 

/ rv n 2 \ 1 /( Q i- Q 2) 
e* = ( . (21) 



The matrix of second derivatives of the effective potential at the critical point (8) is in the 
present case simply a number, and is given by M = — 2aio; 2 . Given that a\ and a<i are 
assumed to have opposite signs, we see that M > 0. 



Perturbative near horizon analysis 

In their study of the attractor mechanism in Einstein gravity, the authors of [5] were able to 
follow two routes towards finding solutions with non-constant scalars. First, they analytically 
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studied solutions that are perturbatively close to double extreme solutions. In this way, 
solutions were found in which the asymptotic values of the scalar fields differ by small 
amounts from their attractor values at the horizon. Reference [5] also constructs solutions 
to the exact equations of motion using numerical techniques. These solutions display the 
attractor behavior over a wide range of asymptotic values for the scalar field. 

In the case of Gauss-Bonnet gravity, the linearized equations for the scalar field in the 
background specified by the metric function (14) cannot be solved in closed form 8 . Therefore, 
to demonstrate the attractor mechanism we will turn very shortly to numerical techniques. 
However, in order to appropriately fix initial conditions for the scalar field near the horizon 
in our numerical work, we first consider the linearized equation for in this region, which 
does yield a simple closed form solution. 

Although for our numerical work below we will specialize to the single scalar model 
described above, this near horizon analysis may be carried through in the general case. Let 
4>i be the critical point values of the scalar fields in a double extreme black hole solution. 
We consider perturbations 

J (r)=0, + e0 /1 (r), (22) 

with e < 1. In order to simplify the analysis, we assume that the fields 0/ are eigenvectors 
of the matrix Mjj with corresponding eigenvalues (5 2 . In the near horizon regime, we have 
b(r) = b H and a(r) 2 = (r — b H ) 2 /R 2 . Setting x = R — b H as in (11), the linearized scalar 
field equations are given by 

3 2 R 2 <b 

x 2 dl<\>i i + 2xd x <j>! i - 2W 1 = • ( 23 ) 
The solutions of this equation are given by 

4>„ = C,(£f (24) 



where of — ( — 1 ± y 1 + 23 f R 2 / b%) /2 and the Cj are arbitrary constants. At this pertur- 
bative level, the attractor mechanism works for a given scalar field 0/, only if the perturbation 
0/i vanishes at the horizon x — 0. We see that this will be the case if the exponent in (24) 
is postive [5]. The exponent erj is always negative and leads to scalar perturbations that 
diverge at x — 0. If (5f > 0, then the exponent 7/ = erf is positive and perturbations with 
this exponent exhibit the attractor mechanism. However, if 3 2 < then no perturbations 
exhibiting the attractor behavior exist. The case 0\ — is discussed in specific class of 
examples in reference [24]. In this case, one must go to higher order in perturbation theory 
to determine the nature of the scalar perturbations, and whether they exhibit the attractor 
behavior. We will assume in the following that all the eigenvalues 3] of the matrix Mjj are 
positive. 

In order to establish initial conditions for all the relevant degrees of freedom, we also 
consider the perturbation equations for the metric functions a 2 and b 2 in the near horizon 

8 Note that corrections to the metric functions would not enter until second order in perturbation theory, 
because of the quadratic nature of the scalar kinetic term and the fact that we are expanding about a critical 
point of V e g. 
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regime. As noted above, these begin at second order in perturbation theory, and hence we 
expand a(x) = a(r) + e 2 a 2 (r) and b(x) = b(r) + e 2 b 2 (r). From equation (10) we can then 
obtain the equation for b 2 (x) 



b 3 H Ei (d^nf 



6R 2 



(25) 



which after insertion of the solution (24) for the first order perturbations to the scalar fields 
can be integrated to give 

C 2 b%(x/b H ) 2 ^ 
^ = 127, (27, -I)**" (26) 

We see that b 2 (x) vanishes at the horizon x — 0, and hence the metric function b approaches 
its attractor value &„ at the horizon. 

Now one can use the third equation in (9) to solve for the purturbation a 2 (x) near the 
horizon and get 



<9 2 \aa 2 + 



b H a 2 



b 2 
°H 



AR 2 7 AR 2 
The integration of this equation gives 



b 7 H 



2b% 



(27) 



aa 2 + 



b H a 2 



AR 2 



x 



-E 
i - 



b\C 2 



r4i? 2 (2 7/ + l)(2 7 , + 2)^ 
b 2 H \ 7? V 



x 



&R 2 'R 2 3^7,(27,-1)^ ' 2b% 



+ 



x 



27/ +2 



(28) 



from which we see that the perturbation to a(r) also vanishes as x — > 0, consistent with the 
attractor phenomenon. 



Numerical solutions 



We now return to the single scalar field model and carry out a numerical analysis of solutions 
to the full nonlinear field equations. From the ij components of the field equations in (9) 
together with (10), we obtain the radial evolution equations 



0"(r) 
b"(r) 
a"(r) 



-2 6 J 







12 



3 b 2 + Aa(l-a 2 b' 2 ) 

1 h_ a V 2 -& 2 a ,2 -2a& (2a'b' + ab") 



a[6 2 + 4a(l-a 2 6' 2 )] 



- Aaa' 2 (l-3a 2 b' 2 )+8aa 3 a'b'b"-b 2 a 2 (d r( p) 2 + 



Kff(0) 



(29) 



We integrate these equations numerically using the Rung-Kutta method. As discussed in [5] , 
one cannot impose boundary conditions near r = oo, because the growing mode of the scalar 
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field in (24) will lead to divergent results near the horizon. Instead, we integrate outward in 
radius. We start the numerical integration at an initial point close to the horizon and use 
the perturbative near horizon results in (24), (26) and (28) to fix initial conditions there. 
We denote the proximity to the horizon by the parameter 



5r 



n - b 



H 



(30) 



With a single scalar field, the strength of the perturbation near the horizon is determined 
by the choice of the single constant C in equation (24). 

As we integrate the system of equations, we check that the numerical solution satisfies 
the constraint given by the rr component of the field equations (9) which contains no second 
order radial derivatives 



•/( r ) =ab' 



a 2 b 2 



(l-a 2 6' 2 ) + (0'r 



Kff(0) 
a 2 6 6 



/3b' 6 2 + 4a(l 



a 2 b' 2 ) 



(31) 



We find that in all cases the constraint is indeed satisfied to appropriate numerical accuracy. 

As a further check, we have also numerically integrated the simpler lowest order per- 
turbative equations for 0i(r), ai(r) and bz(r) expanded about the full, asymptotically flat, 
double extreme black hole background. These perturbative results should be accurate for 
sufficiently small values of the parameter C. We can then check that our results for inte- 
grating the full nonlinear field equation match up with the perturbative results for small 
C. 

Our numerical results are displayed in figure (1). The vertical line in the plots denotes the 
horizon radius bn- The plot of 0(r) clearly shows the attractor behavior. By construction, 
the plots start near the horizon with small deviations determined by the value of C from the 
attractor value 0, which is indicated by the horizontal line on the plot. We see that for large 
radius, the scalar field approaches constant values that can differ quite significantly from 
its attractor value near the horizon. The plot also includes for comparison the results of 
integrating the first order perturbative equation for 0. We see that these results have the 
same qualitative features as the results from the full noninear equations and that the two 
sets of results have good quantitative agreement at the smallest value of C displayed. 

The plot of the metric functions a(r) and b(r) have forms consistent with asymptotic 
flatness, with a(r) going to a constant and b(r) growing linearly at large radius. We will see 
below that their detailed asymptotic behavior gives asymptotically flat solutions with finite 
ADM mass. 



Black hole mass 

In this section we evaluate the ADM mass for our numerical solutions. It is a basic feature 
of the supersymmetric attractor mechanism that the ADM mass is minimized, for fixed 
values of the electromagnetic charges carried by the black hole, when the scalar fields take 
their attractor values throughout the spacetime [3]. It was shown in [5] that this continues 
to be the case for nonsupersymmetric attractors in 4-dimensional Einstein gravity. In the 
following we show that this behaviour persists, at least in our 5-dimensional example, when 
the Gauss-Bonnet interaction term is added to the gravitational Lagrangian. 
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As noted above, our numerical results show that at large radius b(r) increases linearly 
with radius and a(r) approaches a constant value. If we assume that the approach of a(r) 
to its asymptotic value near infinity is power law with some exponent, then for large r we 
have 

Kr) = /r, „>(r) = «i-^, (32) 

where /, and n are constants. For an asymptotically flat spacetime in 5-dimensions, we 
should have n — 2, and Madm in (32) would then be the properly normalized ADM mass. 
Rescaling the time coordinate so that gu approaches —1 at infinity, the asymptotic form of 
the metric is then 

ds 2 = - (l - M \ DM / n ) dt 2 + jJlL + yW 3 . (33) 

where y = jr. It is straightforward to check that the constants / and satisfy the relation 
/doo = 1 in our numerical solutions. In order to check that n = 2 in our solutions and to 
determine the value of Madm, we made a log vs. log plot of — (a 2 — a 2 ^) versus y. Via the 
relation 

log (a%o -a 2 ) = log (^ MA q^/ j - n log (y) , (34) 

we see that for power law falloff of a(r) to its asymptotic value this plot should approach 
a straightline with slope —n. In this manner, we determined that n = 2 in our numerical 
solutions to good accuracy and obtained values of Madm over a range of values of C. 

Figure (4) shows plots of Madm with varying asymptotic values for for a = (i.e. 
pure Einstein gravity) and a — 1. We see that the results are qualitatively similar for the 
two values of a and that in both cases the mass is minimized when the scalar field takes its 
attractor value throughout the spacetime. 



5 Non-extremal black holes 

Finally, we consider nonextremal black hole solutions. We saw in section (3) that nonex- 
tremal solutions exist with the scaler fields fixed at a critical point of V&{4>). The metric 
functions are given by 



a 2 (r) = 1 + 



4a \ 



4a 



2ar 2 



b(r) = r , 



(35) 



where r± are the inner and outer horizon radii with 



r\ = Mil M 2 



v cS (4>io) 



(36) 



We want to ask whether there exist nonextremal attractor solutions in which the scalar fields 
vary between their attractor values at the outer horizon and independent values at infinity? 
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In the Einstein case [5, 18] such solutions do not exist, and we expect to find similar results 
after adding the Gauss-Bonnet interaction. We address this question both analytically, by 
lookiing at perturbations to the scalar field in the near horizon region, and numerically by 
looking at solutions to the full nonlinear field equations. 

We give the perturbative results first. Near the outer horizon the leading order behavior 
of the metric functions expanded in terms of r — r + is given by 



a 2 (r) 



p(r + ,r-)(r — r + ) b(r) ~ r + 



(37) 



where p(r+,r_) = 2(r 2 — r 2 )/(4ar + + r+). The first order perturbative equation for the 
scalar field in the near horizon region is then given by 



(r - r+)<M " + ^ 



P 2 



1r\p 



•Hi 



0, 



(38) 



where ft 2 = didjV e s{4>io)- The solutions for linearized perturbations of the scalar field are 
then given by 



Mr)=CI 




+ DK 




(39) 



where I and K are the modified Bessel's functions of the first and second kind respectively. 
Since K is singular at r = r + we set D — 0. However, we also have Jo (a; = 0) ^ 0. small 
perturbations to the scalar field therefore necessarily modify its value at the horizon, and 
we see that the attractor mechanism no longer holds for nonextremal black hole horizons. 
Figure (3) displays our numerical results for nonextremal solutions to the full nonlinear field 
equations for different values of the parameter C. The fact that the values for at the 
horizon differ from the attractor values agrees with the result of the perturbative analysis 
and shows that the near horizon results extend to nonextremal asymptotically flat solutions. 
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Appendix 



In this appendix, we present further details of the formalism of section (2). The Lagrangian 
is given in equations (1) and (2). The equations of motion for the metric, gauge and moduli 
fields are given by 



1 



Q,u = 2(9 M / )^0 / - g^idMW^ 1 + fM) ^F; x F b v x - -g, v F; v F^' 



(40) 

(41) 
(42) 
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with = + aGj£ B \ where G^ v is the Einstein tensor and G^ 5 * 1 is its Gauss-Bonnet 
counterpart given by 

The nonzero components of the Riemann tensor for the spherically symmetric metric (3) are 
given by 

Rrtr* = -(- + ^t), Rrir^-iT + ^T) 6 ^ ( 44 ) 

a a z b ab 

Rut 3 = A^Wi, Rijk l = {l-a 2 b' 2 ){ liklj l - ljkli l ) (45) 

The nonzero components of the Einstein tensor G^ v and its Gauss-Bonnet counterpart Gj^ B ^ 
are given by 

G tt = ^(l-aW)-Sa^ + *£), G„ = -3-^(1 - a%> 2 ) + 3^, (46) 
Gij = (-(l-a 2 6' 2 ) + 6 2 a 2 (^ + ^) + 2aV(^ + 2^ 7jJ , (47) 

while those of its Gauss-Bonnet counterpart are given by 

C,f B > = -12^(1 -a 2 !/ 2 )(y + ^), G^' = 12^(1 (48) 
or = ^{ (1 - a V^ + ^)-2^<)^ + ^)} 7w . (49) 
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Figure V. This figure displays the results of numerically integrating the field equations with 
C = {5.0, 10, 15}. The parameters of the scalar field model are taken to be a = 1.0, 
Qi = l/y/2 and Q 2 = V% , an = -a 2 = 2.0, which gives b H = 0.904 and = 0.347. The 
numerical integrations are started at 5r = 0.01. In the plot of <ft, the dark lines are the 
result of integrating the full nonlinear field equations, while the lighter lines show the first 
order perturbative results. We see that these results agree well at the smallest value of C 
displayed. The graphs for a and b show that the exact solutions are singularity free and 
asymptotically flat. 
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Figure 2: The black hole mass Madm versus the asymptotic value of the scalar field 0^ for 
a = and a = 1.0. We chose charges Qi = 1/V2 and Q 2 = a/2 , «i = —a 2 = 2.0 and 
5r = 0.01. Plots show that mass increases with 0^. The minimum value of mass is that 
of the double-extreme black hole obtained by setting 0^ equal to its critical value at the 
horizon. 
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Figure 3: The non-extremal black hole . We choose a — 1.0 M — 1.0, charges Q e \ = l/y/2 
and Q e 2 = \/2 , ct\ = —0:2 = 2.0 and 5r = 0.01. Plot shows that the non-extremal solution 
does not exhibit the attractor mechanism. The scalar field is drawn to different values at 
the horizon for different values at infinity. 
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